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Abstract 

Enumeration of all combinatorial types of point configurations and polytopes is a fundamental 
problem in combinatorial geometry. Although many studies have been done, most of them are for 
2-dimensional and non-degenerate cases. 

Finschi and Fukuda (2001) published the first database of oriented matroids including degenerate 
(i.e., non- uniform) ones and of higher ranks. In this paper, we investigate algorithmic ways to classify 
them in terms of realizability, although the underlying decision problem of realizability checking is 
NP-hard. As an application, we determine all possible combinatorial types (including degenerate 
ones) of 3-dimensional configurations of 8 points, 2-dimensional configurations of 9 points and 5- 
dimensional configurations of 9 points. We could also determine all possible combinatorial types of 
5-polytopes with 9 vertices. 

1 Introduction 

Point configurations and convex polytopes play central roles in computational geometry and discrete 
geometry. For many problems, their combinatorial structures or types is often more important than their 
metric structures. The combinatorial type of a point configuration is defined by all possible partitions 
of the points by a hyperplane (the definition given in (12. ip ). and encodes various important information 
such as the convexity, the face lattice of the convex hull and all possible triangulations. One of the 
most significant merits to consider combinatorial types of them is that they are finite for any fixed sizes 
(dimension and number of elements) while there are infinitely many realizations of a fixed combinatorial 
type. This enables us to enumerate those objects and study them through computational experiments. 
For example, Finschi and Fukuda [53] constructed a counterexample for the conjecture by da Silva and 
Fukuda [IB] using their database [3T]. Aichholzer, Aurenhammer and Krasser [2] and Aichholzer and 
Krasser [J] developed a database of point configurations pQ and showed usefulness of the database by 
presenting various applications to computational geometry [3J |3J H] . 
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Despite its merits, enumerating combinatorial types of point configurations is known to be a quite hard 
task. Actually, the recognition problem of combinatorial types of point configurations is polynomially 
equivalent to the Existential Theory of the Reals (ETR), the problem to decide whether a given polynomial 
equalities and inequalities system with integer coefficients has a real solution or not, even for 2-dimensional 
point configurations [37j [46] . Because of this difficulty, enumerations have been done in the following two 
steps. 

The first step is to enumerate a suitable super set of combinatorial types of point configurations which 
can be recognized efficiently. One of the most frequently used structures has been oriented matroids. 
Oriented matroids are characterized by simple axiom systems and many techniques for the enumeration 
have been proposed. Exploiting a canonical representation of oriented matroids and algorithmic advances, 
Finschi and Fukuda 22, 23 enumerated oriented matroids including non-uniform ones (Definition I2.2[) . 
degenerate configurations in the abstract setting, on up to 10 elements of rank 3 and those on up to 
8 elements for every rank. In addition, Finschi, Fukuda and Moriyama enumerated uniform oriented 
matroids (Definition [2~2l in OM(4,9) and OM(5,9) using OMLIB [21], where OM(r, n) denotes the set 
of all rank r oriented matroids on n elements. Aichholzer, Aurenhammer and Krasser [2], and Aichholzcr 
and Krasser [4] enumerated uniform oriented matroids on up to II elements of rank 3. The enumeration 
results are summarized in Table [1] 

In the second step, to obtain all possible combinatorial types of point configurations, we need to 
extract oriented matroids that are acyclic and realizable. Realizable oriented matroids (Definition 12. 6|) 
are oriented matroids that can be represented as vector configurations and acyclic-ness (Definition 12. 3p 
abstracts the condition that a vector configuration can be associated to a point configuration. While 
checking the acyclic-ness is trivial, the realizability problem is polynomially equivalent to ETR (37J US] 
and thus NP-hard. In this paper, we show that the realizability problem can be practically solved for 
small size instances by exploiting sufficient conditions of realizability or those of non-realizability. 

1.1 Brief history of related enumeration 

The enumeration of realizable oriented matroids has a long history. First, Griinbaum 30, 29 enumerated 
all realizable rank 3 oriented matroids on up to 6 elements through the enumeration of hyperplane 
arrangements. Then Canham [13] and Halsey [32] performed enumeration of all realizable rank 3 oriented 
matroids on 7 elements. Goodman and Pollack [27j [28] proved that rank 3 oriented matroids on up to 
8 elements are all realizable. The enumeration of rank 3 uniform realizable oriented matroids on 9 
elements is due to Richter [JT] and Gonzalcz-Sprinberg and Laffaille [35]. The case of rank 4 uniform 
oriented matroids on 8 elements was resolved by Bokowski and Richter-Gebert [TU]. Bokowski, Laffille and 
Richter resolved the case of rank 3 uniform oriented matroids on 10 elements (unpublished). Recently, 
Aichholzer, Aurenhammer and Krasser [2] developed a database of all realizable non-degenerate acyclic 
oriented matroids of rank 3 on 10 elements and then Aichholzer and Krasser [4] uniform ones of rank 3 
on 11 elements. The enumeration results are summarized in Table [2] 

The enumeration of combinatorial types of convex polytopes also has a long history. The combinatorial 
types of convex polytopes are defined by face lattices (See [29, 51 ). All combinatorial types of 3-polytopes 
can be enumerated by using Steinitz' theorem (49] [50]. We can also obtain all combinatorial types of d- 
polytopes with n vertices using Gale diagrams for n < d + 3 [33] [25]. On the other hand, the enumeration 
of combinatorial types of o?-polytopes with (d + 4) vertices and those of 4-polytopes are known to be quite 
difficult [37] 32] . Griinbaum and Sreedharan [3T] listed all combinatorial types of simplical 4-polytopes 
with 8 vertices, and Altshuler, Bokowski and Steinberg [5] those of simplicial 4-polytopes with 9 vertices, 
and then Altshuler and Steinberg [6 those of non-simplicial 4-polytopes with 8 vertices. The enumeration 
results are summarized in Table [3] 

However, there is no database of these objects including degenerate ones or of high dimensions, 
currently. Many problems in combinatorial geometry remain open especially for high dimensional cases 
or degenerate cases, and thus a database of combinatorial types for higher dimensional or degenerate 
ones will be of great importance. For example, characterizing the /-vectors of d-polytopes is a big open 
problem for d > 4 while the same questions for 3-polytopes and for simplicial polytopes have already 
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Table 1: The numbers of simple oriented matroids (n: the number of elements, r: rank) (reorientation 
class, the numbers enclosed by brackets are those of uniform oriented matroids) [H [4j [TO] [13l [22l [23l [26l 
[301 [211 EH [41] 
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Table 2: The numbers of simple realizable oriented matroids (n: the number of elements, r: rank) 
(reorientation class, the numbers enclosed by brackets are those of uniform realizable oriented matroids) 

[H [H HD3 [HI I2S1 E01 12H1 E21 HI] 

been solved [HUHg?]. 

Since Finschi and Fukuda developed a database of oriented matroids [2TJ [22] containing non-uniform 
ones, the realizability classification of larger oriented matroids including non-uniform case has begun. 
Existing non- realizability certificates such as non-Euclideanness [181135) and biquadratic final polynomials 
[5] and existing realizability certificates such as non-isolated elements [33] and solvability sequence [T^] 
were applied to OM(4,8) and OM(3, 9) [531 (Ml HO]- A new realizability certificate using polynomial 
optimization and generalized mutation graphs [40) and new non-realizability certificates non-HK* [24] 
and applying semidefinite programming [36] were proposed and applied to OM(4, 8) and OM(3, 9). Results 
of those classifications are summarized in Figure Q] 

It is important to observe that there are 4803 oriented matroids in OM(4, 8) and 8548 oriented 
matroids in OM(3, 9) whose realizability was previously unknown. 
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Table 3: The numbers of combinatorial types of convex polytopes (the numbers enclosed by brackets are 
those of simplicial polytopes) (n: the number of vertices, d: dimension) [5, 6, 30] 131] 
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Figure 1: Previous classifications of OM(4,8) & OM(3,9) w.r.t. certificates [24 j l39 j l40 l 136] 



1.2 Our contribution 

In this paper, we complete the classification of OM(4,8), OM(3,9) and OM(6,9) (Theorem ITTT|) by 
providing a new method that could successfully find realizations of all previously unclassified oriented 
matroids. 

As mentioned above, the realizability problem is as hard as solving general polynomial inequalities 
asymptotically. There are several methods to solve general polynomial inequalities such as Cylindrical 
Algebraic Decomposition [T5], but the problem size which can be practically dealt with is quite limited, 
and our instances turn out to be intractable. One of the reasons is that those methods compute the 
complete description of cylindrical decomposition of the solution set, which is not necessary for our pur- 
pose. It suffices to find one solution of the polynomial system to decide realizability, which is usually a 
much easier task. It is recently reported that methods based on random realizations are applied success- 
fully to the classification of the realizability of uniform oriented matroids [4] and that of triangulated 
surfaces [33]. However, those methods are not directly applicable to non- uniform oriented matroids. 

In this paper, we take a fresh look at the solvability sequence method [12], which detects the realizabil- 
ity of a given oriented matroid, provided one can eliminate all variables in the polynomial system using 
a certain elimination rule. We extend the elimination rule and introduce some additional techniques so 
that they can be applied to a broader class of oriented matroids. We also use random realizations when 
there are remaining variables in the final step. Using this method, we manage to realize all realizable 
oriented matroids in OM(4, 8), OM(3, 9) and OM(6, 9). This in turn proves that every non-realizable 
oriented matroid in these classes admits a biquadratic final polynomial certificate. 

Theorem 1.1 

(a) Among 181, 472 oriented matroids in OM(4, 8) (reorientation class), 177, 504 oriented matroids are 
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realizable and 3, 968 are non-realizable. 

(b) Among 461, 053 oriented matroids in OM(3, 9) (reorientation class), 460, 779 oriented matroids are 
realizable and 274 are non-realizable. 

(c) Among 508, 321 oriented matroids in OM(6, 9) (reorientation class), 508, 047 oriented matroids are 
realizable and 274 are non-realizable. 

We note here that the classification of OM(6, 9) is obtained from the classification of OM(3, 9) and the 
duality of oriented matroids [8]. As a byproduct, we obtain the following results. 

Theorem 1.2 

(a) There are exactly 15, 287,993 combinatorial types of 2-dimensional configurations of 9 points, 
105, 128, 749 5-dimensional configurations of 9 points and 10, 559, 305 3-dimensional configurations 
of 8 points. 

(b) There are exactly 47, 923 combinatorial types of 5-dimensional polytopes with 9 vertices. Among 
them, 322 are simplicial and 126 are simplicial neighborly. 

Our classification results with certificates are available at 

http : //www-imai . is . s . u-tokyo . ac . jp/~hmiyata/oriented_matroids/ 

To make the results as reliable as possible, we recomputed realizability or non-realizability even for ori- 
ented matroids whose realizability had been known previously. In the above web page, realizations of all 
realizable oriented matroids and final polynomials of all non-realizable oriented matroids are uploaded. 
One can check correctness of our results there. 

Organization of the paper: 

In Section 2, we present some basic notions on oriented matroids. Then we discuss a standard method 
to find realizations in Section 3. We first review the existing method to decrease the size of a polyno- 
mial system, and explain a new method to search for solutions of polynomial systems. We apply these 
methods to the classification of OM(4,8), OM(3,9) and OM(6,9) in Section 4, and conclude the paper 
in Section 5. 

2 Preliminaries 

In this section, we review basic notions on oriented matroids that are used in the paper. For further 
details about oriented matroids, see [8]. 

2.1 Point configurations and their combinatorial abstractions 

Let P = {pi, . . . ,p n } be a point configuration in R r_1 . We define a map x '■ {!>•••) n Y ~^ {+? — ■> 0} by 

X(h, ■ ■ ■ , v) ■= sign(det(u il , . . . , v ir )), 

where v\ :— {^-^j ■>■ ■ ■ > v n '■= e ^ are ^ e associa ted vectors of pi, . . . ,p n . We define the map \ as 

the combinatorial type of P, which satisfies the following properties (a), (b) and (c) with E = {1, . . . , n}. 

Definition 2.1 (Chirotope axioms) 

Let E be a finite set and r > 1 an integer. A chirotope of rank r on E is a mapping \ '■ E r — > {+1, —1, 0} 
which satisfies the following properties for any i\, . . . , i r , j\, . . . , j r G E. 

(a) x is n °t identically zero. 
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(b) x(v(i) , • • • , i a {r)) = s g n ( <T )x(*ii • • • > V) f° r a ll Hi • • • > V € E and any permutation cr. 

(c) If X(js,i2,---,i r ) ■ x(ji, ■ ■ ■ , 3s-i,h,js+i, ■ ■ ■ , 3r) > for all s = 1, ...,r, then 

X(«l,---,«r) •■•,>) > 0. 

We note here the third property is an abstraction of Grassmann-Pliicker relations: 

r 

[il... i r ]\jx . . .j r ] - ^ \jsh ■ ■ ■ ir]\jl ■ ■■js-liljs+l ■ --jr] = 0, 
s=l 

where . . . ,i r ] := det(w il , . . . , for all ii,...,i r G E. We define an oriented matroid as a pair 
of a finite set E and a chirotope \ '■ E r — >• {+1,-1,0} satisfying the above axioms. From now on, 
we set E := {l,...,n} throughout this section. Since \ is completely determined by the values on 
A(n, r) := ...,v) G N r | 1 < ii < • • • < i r < n}, we sometimes regard the restriction xU(n,r) as the 
chirotope itself. We call the pair (E, {\, — x}) a rank r oriented matroid with n(= \E\) elements. The 
concept of degeneracy is also defined for oriented matroids as follows. 

Definition 2.2 An oriented matroid (E, {x, — x}) is said to be uniform if ■ ■ ■ i^r) ^ for all 
(ii, . . . , i r ) £ A(n, r), otherwise non-uniform. 

Note that every (d + l)-subset of a (i-dimensional point configuration spans a o?-dimensional space if and 
only if the underlying oriented matroid is uniform. 

For a point configuration P = {pi, . . . ,p n } C M. d , the data of chirotope is known to be equivalent to 
the following data. 

V£ := {(sign(a T pi - b), sign(a T p 2 - 6), ... , sign(a T p„ - 6)) | a e R d , b G M}. (2.1) 

An element of Vp is called a covector. Another axiom system of oriented matroids ( Covector axiom) can 
be obtained by abstracting properties of V* P . See [H Chapter 3], for details. 

Note that every point configuration has a covector (+ + •••+). Abstracting this property, acyclic 
oriented matroids are defined. 

Definition 2.3 An oriented matroid is said to be acyclic if it has a covector (+ + ••• +). 

It is known that there is the one-to-one correspondence between acyclic realizable oriented matroids and 
combinatorial types of point configurations. 

From Vp, we can read off the convexity of P. For i = 1, . . . , n, Pi is an extreme point of P if and 
only if there is a covector (+••• + +••• +) G Vp. In this way, matroid polytopes are introduced as 

i-th component 

abstractions of convex point configurations. 

Definition 2.4 An acyclic oriented matroid on a ground set {1, . . . , n} is called a matroid polytope if it 
has a covector (+••• + + •••+) for all i = 1, . . . , n. 

z-th component 

For a matroid poltyope, its facets are defined by non- negative cocircuits i.e., non-negative minimal cov- 
ectors. For details on matroid polytopes, see [8] Chapter 9]. 

2.2 The realizability problem 

Every vector configuration has the underlying oriented matroid, but the converse is not true because 
"non-realizable" oriented matroids exist. 
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Definition 2.5 (The realizability problem of oriented matroids) 

Given a rank r oriented matroid M = (E, {x, —x}) with n elements, the realizability 'problem for 
(E, {x, — x}) is t° decide whether the following polynomial system has a real solution V\, . . . , v n £ M. r : 

sign(det(wj l; . . . ,v ir )) = x(h, ■ ■ • , v) f° r all (h, ...,i r ) £ A(n,r). (2.2) 

Definition 2.6 A rank r oriented matroid is said to be realizable if it arises from an r-dimensional vector 
configuration, otherwise non-realizable. 

Not every realizable oriented matroid admits a point configuration because a positive combination of 
some vectors can be while the same is not true for the associated vector configuration of any point 
configuration. We can extract combinatorial types of point configurations by picking up acyclic realizable 
oriented matroids. 

2.3 Isomorphic classes 

In this paper, we consider only simple oriented matroids, those without parallel elements and loops, 
see [5]. For simple oriented matroids, we consider the following two equivalent classes. 

Definition 2.7 Let M = (E, {x, -%}) and M' = (E, {%', -%'}) be oriented matroids. 

(a) M and M' are relabeling equivalent if 

X(k, . ■ • , v) = x'(4>(n), ■ • ■ , <Mv)) for all i%, . . . ,i r £ E 

or 

x(ii, . . ■ , v) = -x'(4>(n), <f>(i r )) for all i%, . . . ,i r £ E 
for some permutation on {1, . . . , n}. 

(b) M and M' are reorientation equivalent if M and —aM' are relabeling equivalent for some A Q\ E, 
where —aM' is the oriented matroid determined by the chirotope — ax' defined as follows. 

-Ax'ih, ...,i r )~ (-l) |An{il '-^ }l x'(H, • • • , ir) for h, . . . , i r £ E. 

We note here that the realizability is completely determined by the reorientation classes of oriented ma- 
troids, the equivalence classes defined by the reorientation equivalence. This is because if any oriented 
matroid in a reorientation class is realizable, then every oriented matroid in the class is realizable. A 
database of oriented matroids by Finschi and Fukuda [21] consists of the representatives of the reorienta- 
tion classes. On the other hand, we say that point configurations P, P' has the same combinatorial type 
if the combinatorial type of P and that of P' belong to the same relabeling class. 

3 Methods to recognize realizability of oriented matroids 

Recognizing that a given oriented matroid is realizable amounts to finding a solution of the associated 
polynomial system (|2.2[) in Definition 12.51 

Our strategy is as follows. We first simplify the polynomial system as much as possible, namely, 
eliminate as many variables as possible using simple elimination rules and then try random realizations 
if no further elimination of variables is possible by the elimination rules. 

3.1 Inequality reduction techniques 

There are three critical parameters to measure difficulty of solving polynomial systems: the number of 
variables, the degrees of variables and the number of constraints. We must try to keep each of them low. 
For such purposes, we employ some techniques used in [Til H21 EH HQ]. Let us review these techniques 
briefly. 
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3.1.1 Fixing a basis 



A technique explained in the following was introduced in [T^] and was used to reduce the degrees of 
constraints and the number of variables in [T^l IM1 EH] • 

We assume that x(ii,...,i r ) = + for some ii,...,i r by taking the negative of x if necessary. Let 
My '.— (vi,...,v n ) G R dXTl be the representation matrix of a vector configuration V. Because the 
combinatorial type of V is invariant under any invertible linear transformations, we can assume that 
v bl = (1,0, ...,0) T ,v b2 = (0, 1,0, ...,0) T , ...,v br = (0, ...,0, 1) T for an r-tuple (&i, ...,b r ) G A(n,r) such that 
x{bi, ...,b r ) = +. We call such an r-tuple of indices a basis. We obtain a new polynomial system as 
follows. 

sign(det («,!,..., w» r )) = x(h,—,i r ) for all (ii, v) G A(n,r), and 
u 6l = (1,0,...,0) T ,^ 2 - (0,l,0,...,0) T ,... jt ; 6r = (0,...,0,1) T . 

The resulting polynomial system (13. ip depends on the choice of bases. In the next section, we present 
how to find a suitable basis. 

Finally, note that the degree of each constraint can be computed easily by the following formula: 

deg(det(v il , ...,v ir ))) = ....,b r } \ 

and that the sign of each variable Vki is determined by the obvious equation: 

v k i = det (v bl , . . . , Vb k _ t , vi , Vb k+1 , Vb r ) ■ 



3.1.2 Minimal reduced systems 

Most of the techniques in this section are introduced in [12] for uniform oriented matroids to reduce the 
number of constraints and are extended to general oriented matroids in |38[ 140] . 

Recall that in the polynomial system (|3.1j) . which describes the condition of realiz ability, there are 
("J constraints. There are some possible redundancies in these constraints, as one may reconstruct x 
according to partial values of x by using Axiom (c) of a chirotope. For example, if x(l, 2, 3) = x(l, 4, 5) = 
x(l, 2, 4) = x(l, 3, 5) = x(l, 2, 5) = +, then wc obtain x(L 3,4) = + using Axiom (c). For a subset R of 
A(n, r), we denote (R) the set of all r-tuples whose x signs are implied by the sign information on x over 
r-tuples in R and the chirotope axioms. A minimal subset R(x) °f A(n,r) needed to decide x is called 
a minimal reduced system for (E,{x,~x})- First, we observe that generalized mutations to be defined 
below must be contained in every minimal reduced system of x- 

Definition 3.1 An r-tuple A G A(n, r) is called a generalized mutation of x if there exists an oriented 
matroid (E, {x', — x'}) such that x(^) — x'(m) f° r an M h-{ n i r ) \ an( i xW x'W- 

This definition is different from that of [40l|38], but turns out to be more natural because of the proposition 
below. 

Let us denote the set of all generalized mutations of x by GMut(x)- A nice characterization of the 
mutations of uniform oriented matroids is given in [451 Theorem 3.4.], but such a characterization is not 
known for generalized mutations. However, one can compute GMut(x) using the following proposition, 
which is an immediate extension of |45i Proposition 3.3.]. 

Proposition 3.2 An r-tuple A := (ix, . . . ,i r ) £ A(n, r) is a generalized mutation of x if and only if A is 
not determined by Grassmann-Pliicker relations, i.e., the following condition holds: 

X(h,---,3r) =0, or 

X(h,-- ■ ,jr) + and {x(j»> i 2, ■ ■ ■ , i r )x(h> ■ ■ ■ > js-i, h,ja+u • • • »Jr) I s = l,...,r} D {+,-} 
for all 1 < ji < ■ ■ ■ < j r < n. 
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Proof. In this proof, we use the notation (fci, . . . , fc r ) to represent an r-tuple (k[, . . . , k' r ) such that 
{k[, . . . , k' r } = {fci, . . . , k r } and k[ < ■ ■ ■ < k' r , for fci, . . . , k r G N. 

To prove the if-part, let us take a map x : A(n, r) — )• {+,—,0} satisfying Condition p. 21) . We 
consider a map x' '■ {!>..., n} {+, ~ , 0} satisfying Axioms (a) and (b) of a chirotope and x'(m) = 
x(/«) for all fi e A(n,r) \ {A} and xW x'W- We prove that x' also satisfies Axiom (c). For 
(fci, ... , fc r ), (Zi, . . . , Z r ) € A(n, r), assume that 

x'(Z s , fc 2 , ■ • • , fc r ) • x'(h, ■ ■ ■ ,Z s _i, fci,Z s+ i, . . . ,Z r ) > for all s = 1,. ..,r. (3.3) 

Under this condition, we prove x'(fci5 • ■ • ? k r )x'(h, ■ ■ ■ , Z r ) > by the following case analyses. 

(I) fci e {Z!,...,Z r }. 

Let t 6 {1, . . . , r} be an integer such that k\ = It- 

X(h, k 2 , ■ ■ ■ , K) ■ x'(h, lt-i,h,l t +i, • • • , Z r ) = x'(fci> • • • , K) • x'(h, . . . , Z r ) > 0. 

(II) fci g {Zi,...,Z r }. 

(II-A) |{fc 1 ,...,fc r }n{Z 1 ,...,Z r }| =r-l. 
Let us take Z u ^ {fci, . . . , fc r }. Then we have 

X'(lu,h, • ■ • , fcr) ' X'('l, • • • > fci, Wl, • • • , Zr) = (-1)" • x'Gl, • • • , It) ' (-1)" • x'(fcl, • ■ • , fcr) 

= X'(Zl,...,Zr)-x'(fcl,...,fcr) >0. 

(II-B) |{fc 1 ,...,fc r }n{Zi,...,Zr}| <r-l. 

(i) {fc x ,...,fcr} = A. 

Since fci {Zi, . . . ,l r } and |{fci, . . . , fc r }n{Zi, . . . ,l r }\ < r-1, (Z SJ fc 2 , . . . , fc r ) ^ A and (Z x , . . . ,Z s _i, fci,Z s+ i,. . . ,Z r ) 
A for s = 1, . . . , r. Therefore, Condition (|3.3[) implies 

X(Z S , fc 2 , ■ ■ ■ ,fc r ) • x(h, ■ ■ .,l a -uh,ls+l, ■•■,l r ) > for all s = 1, . . . ,r. 
To satisfy Condition (|3.2[) . the following must hold. 

X'(Zl,...,Zr)=x(Zl,---Jr)=0. 

(ii) (h,...,l r ) = A. 
Proved similarly to Case (i). 

(iii) (fci, . . . , fc r ) ^ A, (Zi, ...,l r )^\. 

First, we consider the case when there exists so such that A = (Z SQ , fc 2 , . . . , fc r ) and x(h, • • • i ^s -i5 fclj ^s +i> • • • , lr) 7^ 
0. Let us write down Condition (|3.2[) for x under ii := l So , i 2 := fc 2 , . . . , i r :=k r and ji := Zi, . . . ,j So -i '■— 

lso — lj jso ■ — fci j Jso + 1 ■ — fcso + 1 5***5 3t ' lr ■ 

{X(l s ,k 2 , ... ,k r )x(h, ■■■ Jso-l^iJso + l, ■■■ ,l s -l,l So ,l s +l, ■■■ Jr) | S = 1, ...,8q - 1} 

U {x('s5 fc2,** ■ ,k r )x(h, ■ ■ . ,Z s _i,Z SQ ,Z s+ i, . . . ,Z So _i,fci,Z So+ i, . . .,Z r ) I s = s + 1, ...,r} 

U {x(fcl,fc 2 , ■ • ■ ,fcr)x(jl, ■ • ■ ,jr) I S = l,...,r} 

{(-1) • X(ls,h, • • • , fcr)x0l5 • ■ -,^0-1,^0,^0 + 1 • --ls-1, fcl,Z s +l, . ..,Z r ) | S = 1, S - 1} 

U {(-1) • x(ls,k 2 , . . . , fcr)x('i5 • • • ,Z s -i,fci,Z s+ i . . . Z So _i,Z So ,Z So+ i, . ..,Z r ) I s = so + 1, 

U{x'(fcl,...,fcr)x'(Zl,...,Zr)} 

^{ + 5-} 

This implies x'(fci, • ■ • , k r )x'(h, ■ ■ ■ ,l r ) = +■ 

We consider the other case. In this case, Condition Q3.3P implies 

X(Z S , fc 2 , . . . , fc r ) • x(<4> ■ • • , fci, ^+1, ■ • • , Zr) > for all s = 1, . . . , r 
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and thus 

X'(fci, . . . , k r )x'(h, ...,l r )= x(fci, . . . , k r )x{h,. ..,l r ) = +. 

Finally, we conclude x'(^i, • • • , k r )x'(h> ■ ■ ■ Jr) > for all cases. This proves the if-part. 

The only if-part is proved by contradiction. Suppose that there exist 1 < j% < • • • < j r < n such that 

xtii, ■ ■ ■ , jr) ^ and {x{j s , «2, • • • ,V-)x(ji, ■ • ■ h,j s +i, ■ ■ ■ Jr) | s = 1, ...,r} = {+} (or {-} or {0}). 

In this case, we obtain x(h, ■ ■ ■ , v) = x(ji, ■ ■ ■ , 3r)x(Ju *2, • • • , ir)x{h, h, ■ ■ ■ ,jr) using Axiom (c) of a 
chirotope. On the other hand, for a map x' : A(n,r) — > {+,—,0} taking the same value as that of x 
except for A, the following holds. 

X'(il, ...,i r )^ X{h, ...,ir) = X(jl,- ■ ■ ,jr)x(h,i2, ■ ■ ■ , ir)x(H , h , ■ ■ ■ ,jr) 
= X'til, ■ ■ ■ ,jr)x'(jlj2, ■ ■ .,i r )x'(h,j2, ■ ■ ■ , jr)- 

This means that x' violates Axiom (c) of a chirotope. I 

The above proposition shows that GMut(x) C R(x)- Therefore, we can compute R(x) by the following 
procedure, which is a slightly modified version of the algorithm in [TT1 138] , 

Algorithm 3.3 (Computing a minimal reduced system for x) 

Input: A chirotope x '■ {1) n Y ~> {+, 0}, b E A(n, r) s.t. x(b) ^ 0. 

Output: A minimal reduced system R(x)- 

Step 1: Compute GMut(x) and set R := GMut(x). 

Step 2: C := (R). If C = A{n, r), return R. Otherwise go to Step 3. 

Step 3: Choose /i £ A(n, r)\C that minimizes |/x \ 6|, and add it to R. Go to Step 2. 

As a result, we obtain a new polynomial system as follows. 

sign(det(t!i 1 ,...,u ir )) = x(h,—,ir) for (»i,...,i r ) € R(x)- (3-4) 

One might be able to simplify the reduced polynomial system by selecting a different basis. One way is 
to search a better basis 6', in the sense of the totality of the degrees of the constraints, by computing 
X}/3e_R( x ) 1^' Various weight functions of constraints are considered in (TTJ [38] . This is the subject of 
the next section. 

3.1.3 Normalization (eliminating homogeneity) 



A technique to be explained here is used in [381 1" to reduce the number of variables and the degrees of 
polynomials, and in addition, to eliminate some of the equality constraints. 

First, we negate negative variables to make all variable non-negative, and obtain the following new 
polynomial system: 



sign(det(?4,...,^ r )) = (-l) s *i-*r ■ x (ii, ...,i r ) for (i u ...,i r ) G R(x)- 
v' bi = (1,0,. ..,o) T X 2 =(0,1,0, ...,0) T ,.. .,v' br =(0,...,0,1) T , 



where s^...^ denotes the number of negative variables in vu ly u r iu Vu 2 , ...,v r i r . 

Let ai, a m Pit Pr be arbitrary positive numbers and W\, ...,w r £ R™ be row vectors of My := 
(v[,...,v' n ). Then 

det(oti 1 v' h ,...,a ir v' ir ) = a it ■ ■ -a ir det(v' h , u-J, 



/ PiWi(ii, ...,i r ) N 
det : 

\p r w r (ii, ...,l r )j 



fwi(ii,...,i r 
= Pi ■ ■ ■ p r det : 

\w r (ii, ...,i r )j 



10 



where Wj(ii, i r ) denotes the vector whose fc-th element is an i^-th element of Wj for 1 < i\ < ... < i r < n 
and 1 < j < r. This allows us to fix two indices 1 < i < n and 1 < j < r and to assume that every 
component of the vectors v[ and Wj is or 1. 

Furthermore, we can eliminate some of equality constraints as follows. Let us classify the constraints 
arising from ••-> ir) according to the values of \{b\, b r } \ i r }\. 

Case (a) |{6 1} b r } \ i r }\ = 1. 

The constraint corresponds to the sign constraint "u^ > 0" or "u^ — 0" for certain 1 < i < n, 
l<J<r. 

Case (b) \{b 1 ,...,b r }\{i 1 ,...,i r }\=2. 

Let br} \ {ix, i r } =: {«,«'} and {ii, ...,i r } \ {b±, ...,b r } =■ {j,j'}- Then the constraint 

is of the form v' i4 v'-,-, — v'-,,v'--, — 0, where each variable may be fixed to or 1. If the row i is 
normalized, the constraint become v'-, { , — w', = and we can eliminate the constraint in the trivial 
way without increasing the degree of the polynomial system. If the column j is normalized, the 
constraint become v'^, — v' v y — and can be eliminated in the trivial way. 

Case (c) |{6i,...,& r }\{i 1 ,...,v}| > 2. 

The constraint cannot be eliminated trivially. 

Remark 3.4 Equality constraints are main obstacles for random realizations. The normalization tech- 
nique turns out very useful in removing equations so that random realizations can be applied. In addition, 
the technique decreases the number of variables and that of constraints by 1. Therefore, we minimize the 
number of equality constraints over all choices of bases and columns and rows to be normalized, which 
can be computed easily. 

3.2 Searching for solutions of polynomial systems 

Now we present a practical method to find a solution of a polynomial system. 

Cylindrical Algebraic Decomposition method [14] eliminates variables preserving the feasibility until 
polynomial systems contain only one variable, solves the one-variable polynomial systems, and then lifts 
the projected solutions. It may sound quite simple but each step is a very hard task in general. Therefore, 
we take a fresh look at the solvability sequence method [12] . by which some variables might be eliminated 
in a simpler way. 

Proposition 3.5 (12 ) Let hjhih > be integers and R^' ,Flf ,Lj ,lf?\Pk be polynomials for 
i = = 1, ...,Za, fc = 1, ...,{3. Then the feasibility of rational polynomial system: 

Pk(x u ...,x n ) > (or =) (k = l,—,h) 
is equivalent to that of the following polynomial system: 

\Lf\xi, ...,x n )Rf \xi, ...,i n ) < Rf'(x 1 ,,..,x n )Lf\xx,...,x n ) (i = l,...,h,j = l,...,/ 2 ), fi x 
[P k ( Xl ,...,x n ) > (or =) (k = l,...,l 3 ) 

(2) (2) 

under the condition R\ (x\, x n )Lj (x\, x n ) > for i = 1, li,j = 1, li- 

Note that a solution (j/*, x* t £*) of the original system (|3 . 5[) can be constructed from a solution 
(x*, x* b ) of the resulting system by 

, _ min{R l (x* 1 ,...,x* l ) \ i = l, h} + max{Lj(xl , ...,<) | j = 1,...,Z 2 } 
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This elimination rule is used in the solvability sequence method under the bipartiteness condition for 
determinant systems, i.e., polynomial systems whose constraints are of the form "sign(det(wi 1 , ...,w ir )) = 

(2) (2) 

X(iii •••> if)-" I n the determinant system, one can detect the signs of Rl (xi, ...,x n ) and IA {x\, x n ) 
in advance using the information of x, and can solve the polynomial system by y. 

Definition 3.6 ([12]) 

Consider the polynomial system (|3.5[) arising from a determinant system. Then each constraint can be 
rewritten as "sign(det(«i 1 , Vi r )) — iV)" where (ii,...,i r ) € A(n, r) and can be indexed by 

(xi x ) 

{«!,..., i r }- Let A be the set of indices which defines constraints of the form y < —frr — '"' and B the 

index set whose elements define constraints of the form y > — |^ . The polynomial system (|3.5[) 

is said to be bipartite if ■ ■ ■ , i r } H ■ ■ ■ ,jr}\ = r — 1 for all {i\, . . . , i r } 6 A, {ji, . . . , j r } G 

In [T2] , it is proved that the feasibility of the polynomial system (|3.5p is equivalent to that of the following 
polynomial system under the bipartiteness condition: 

P k (xi,...,x n ) > (or =) (k= l,...,h)- 

Therefore, the elimination does not produce new constraints under the bipartiteness condition, and one 
can proceed with eliminations without creating inconsistency. It detects realizability if all variables are 
eliminated by this elimination rule. However, the restriction of the bipartiteness condition is very strong. 
We can easily get rid of the restriction of bipartiteness condition by allowing the elimination rules to 
destroy the determinant system. This relaxation turns out to be very useful when it is used together with 
branching rules, which will be explained later. 

Before explaining branching rules, we also consider an elimination rule for polynomial systems con- 
taining equalities. 

Proposition 3.7 Let l\, £2, h > be integers and Pi, Ej, E, Qk be rational polynomials for i = 1, — 
1, I2, k = 1, I3. Then the feasibility of rational polynomial system: 

' Pi(x\,...,x n ,y) > (i = 1, ...,/].), 

E j [x 1 ,...,x n ,y) = (j = l,...,h), 

y = E(xi, ...,x n ), 

,Qk(xi,-.,x n ) > (or =) (k = 1, ...,Z 3 ) 
is equivalent to that of the following rational polynomial system: 

Pi(xi, ...,x n ,E(xi, ...,x n )) > (i= 

E j (x 1 ,...,x n ,E(xx,...,x n )) = (j = 1,...,Z 2 ), (3.8) 
^Qk(xi, -,x n ) > (or =) (fc = 1,...,Z 3 ) 

Proof. A solution (x* , x* , y*) of the original system (|3.7[) can be constructed from a solution 
(ac*, x* n ) of the resulting system (13.81) by y* := E(x*, 1*). I 

In Propositions 13.51 and 13. 71 we regard a variable y appearing in these forms as redundant. To apply 
these elimination rules to as many variables of this type as possible, we consider the following additional 
rules, which are called branching rules. 

Proposition 3.8 Let £1,^2 > 0. The polynomial system: 

(Ai(xi,...,x n )y < Bi(xi,...,x n ) (i = l,...,h), 
\Pj(x!, ...,x n ) > (or =) (j = l,...,l 2 ) 



(3.7) 



12 



is feasible if and only if one of the following rational polynomial systems is feasible 

'sign(Ai(a;i, x n )) = s(i) (i = 1, h), 



Bi(x t . 

y Aiixi, 

Bi(xi, 



.) 



y > Ai(x!,...,x n ) 

Bi(xi,...,x n ) > 
Pj(x u ...,x n ) > (or 



=)0 



(i < = +), 

(i < ii,s(i) = -), 
(i < «(*) = 0), 
(7=1,..., la) 



(3.10) 



for 



0}. 



Proof. We partition R n+1 into 3' 1 parts R s := {(a^i, x n , y) | sign(Aj(a;i, x n )) — s(i) (i = 1, ?x)} 
for s : {1, /} — > {+, — , 0}. Let S be the solution space of the polynomial system (13.91) . Then S D R s 
is represented by the polynomial system (|3.10[) . S is non-empty if and only if S n i? s is non-empty for 
some s : {1, 1} —>•{+,—, 0}. I 

Observe that one can actually reduce the search range as follows 



sign(Ai(xi, ...,£„)) = s(i) (i = 1, 

^=1 (i </i,*(») = +), 



y > Aj(xi,...,x„) 

Pj(a;i,...,a; n ) > (or =) (j 



(i < l\, s(i) = 
l,...,/ 2 ) 



(3.11) 



for s : {l,...,/i} -> {+,-}. 

This is because (cc* + ei, . . . , a;* + e„,y* + e n +i) is a solution of (|3.9p for a solution (i*, . . . , a;* , y*) 
of (|3.9p and any sufficiently small ei, . . . , e n > 0. The following is a branching rule for the case when y 
appears in equations. 

Proposition 3.9 Let li,h,h > 0. The polynomial system: 

' Ai(xi,...,x n ,y) > (i = l,...,h), 

Aj(xi, ...,x n ,y) = (J = h + 1, ...,h + k), 
A{x\, ...,x n )y = B(xi, ...,x n ), 

,Pfc(xi, ...,x„) > (or =) (fc=l,...,Z 3 ) 

is feasible if and only if one of the following rational polynomial systems is feasible 

A(xi,...,x n ) = 0, 
B(xx, ...,£„) = 0, 

Ai(xx,...,x n ,y) > (i = l,...,Zi), 

(j 



(3.12) 



^ (xi, x n , y) = 



A(xi,. 
Ai(xi, 
Aj(xi, 
Pk(xi, 

A(x u . 
A i (x 1 , 
Aj(xt, 
Pk{xi, 



..,x n ) > (or =)0 (k 

;X n ) > 0, 



ll + l,...,h + l 2 ), 
l,...,l 3 ) 



B\xi " 
Ajxl 

..,x n ) > (or =) 



(j 
(k 



■ ■''"■■> A(xll..]x") ) — U 



.,X, 

..,x. 

.., X 



< 0, 



> 



»' A(£Cl,...,X„)' 

% :j )=o a 



.,x„) > (or =) 



(i 
(j 
(k 



l,...,h), 

: h + l,...,h+l 2 ), 

l,...,h), 
: II + l,...,h + l 2 ), 

: l,...,l 3 ) 



(3.13) 
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Proof. This is proved similarly to Proposition 13.81 We consider the following partition of R n+1 : 
R s := {(xi,...,x n ,y) | sign(A(x 1 , x n )) = s} for s = +,-,0. I 

To solve the polynomial system, the following operations are applied repeatedly. We first choose a 
variable y that can be eliminated by the above 4 rules. Then the branching rule in Proposition 13.81 or 
Proposition ^. 91 is applied to obtain a set of polynomial system p.lO|) or (|3.13p . We choose a sign pattern 
and apply the elimination rule in Proposition [33] or Proposition l3.7l and check the feasibility recursively. 
If feasibility of the polynomial system for this sign pattern is certified, the original polynomial system 
is proved to be feasible. Otherwise, we backtrack and try another sign pattern until feasibility is certi- 
fied. If no feasible polynomial system is found, we give up deciding feasibility. We adopt the following 
termination condition. If all variables are eliminated and the system is consistent, the original system 
is feasible. If a system does not have variables that can be eliminated, we try random realizations to 
prove the feasibility. There are many polynomial systems, which are trivially feasible but are hard to 
solve algebraically For example, the following polynomial system is clearly feasible but is algebraically 
complicated. 

^100 _ y 49 + 10(% 23 <1)X * > 2j / > 2 . 

Random realizations work well as long as the solution set is sufficiently large and are not affected so much 
by the algebraic complexity. 

Algorithm 3.10 Sol(P) (P: polynomial system) 

1. If there are no variables that can be eliminated by the above 4 rules in P, try random assignments 
to the remaining variables. If a solution is found, return 'feasible.' Otherwise return 'unknown.' 

2. Choose a variable y to eliminate in P. Apply one of the branching rules (Proposition l3.8[ Proposition 
I3.9|) to obtain a set of polynomial systems P{, . . . , P' m . 

3. For i — 1, . . . ,m, apply one of the elimination rules (Proposition 13.51 Proposition I3.7P to P[ and 
obtain a new polynomial system Qi. If Sol(Qj) returns 'feasible,' return 'feasible.' 

4. Return 'unknown.' 

We sometimes generate too many branchings and thus apply the iterative lengthy search to the 
following tree search problem. The root node consists of the original polynomial system. Starting from 
the root node, we expand nodes using the elimination rules and the branching rules repeatedly. We decide 
whether we arrive at goal nodes or not using random assignments to the remaining variables. In this 
setting, we define the cost of each node x by c(x) :— log2{rib), where n& is the number of branching at x, 
and apply the iterative lengthening search to it by increasing the limit of the total cost by 1. 

Remark 3.11 Equalities containing no square-free variables cannot be eliminated by the above rules. 
The probability of yielding a solution to equality constraints by random realizations is and if we cannot 
eliminate all equality constraints, it is highly unlikely that the above method generate a realization. In 
this case, which is quite rare for small instances, we need to use general algorithms such as Cylindrical 
Algebraic Decomposition [14] to find solutions. 

We point out another tractable case, namely, the case when the ideal generated by the equality 
constraints is zero-dimensional (as an ideal in C[x s | s G S], which denote the polynomial ring over C 
in variables x s which appear in equality constraints). In other words, it is the case when the equality 
constraints have a finite number of complex solutions. In this case, we extract real solutions, substitute 
each solution to the original polynomial system and apply the above methods. We can check the zero- 
dimensionality and solve the equalities using Grobner basis [T3] . 
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4 Realizability classification 



We apply our method to OM(4, 8) and OM(3,9). All computations are made on a cluster of 4 servers, 
with each node having total 16-core CPUs (each core running at 2.2 GHz) and 128 GB RAM (each 
process uses only 1 CPU). 

First, we apply the polynomial reduction techniques described in section 3.1 to OM(3, 9). Table [4] 
and Table [5] show the distributions of the number of variables and constraints of the resulting polynomial 
systems for OM(3,9). 



No. of variables 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11-18 


No. of OMs 





21 


1,411 


13,261 


47,888 


91,855 


121,977 


107,869 


59,284 


16,814 


673 



Table 4: Number of the variables (OM(3,9)) 



No. of variables 


1-2 


3-4 


5-6 


7-8 


9-10 


11-12 


13-14 


15-16 


17-18 


No. of OMs 


659 


22,340 


107,465 


168,995 


114,595 


38,944 


7,330 


694 


31 



Table 5: Number of the constraints (OM(3,9)) 

We note here that the case of more than 10 variables occurs as an exceptional case, where the 
polynomial system consists of the following type of constraints. 

xy > zw, x, y, z,w > 0. 

Nakayama |38) proved that oriented matroids with such polynomial systems which do not admit bi- 
quadratic final polynomials are realizable. We detect polynomial systems of this type and stop the 
polynomial system reductions because our method of realizations can solve such polynomial systems 
easily. 

We apply our method to search for realizations to the resulting polynomial systems and manage to 
find realizations of 460,778 oriented matroids in OM(3, 9). Table [6] shows the distribution of the time 
which our method consumed to find realizations. 



Time (ms) 


1-10 


10-10 a 


10 2 -10 3 


10 a -10 4 


10 4 -10 b 


10 b -10 b 


10 b -10 y 


lO'-lO* 


10 8 -10 y 


No. of OMs 





241,593 


163,927 


25,893 


19,559 


7, 534 


1,599 


615 


58 



Table 6: Computation time (OM(3,9)) 

Remaining one oriented matroid turns out to be an irrational one, which was found by Perles (See [30|, 
p. 73]). We can realize it using Grobner basis. 

As a result, we obtain complete classification of OM(3,9). In addition, it leads classification of 
OM(6, 9) because the realizability of oriented matroids are preserved by the duality |S]. Similarly, we 
apply our method to OM(4, 8) and manage to give realizations to all realizable oriented matroids in 
OM(4, 8) (Theorem II. 1| except for two irrational ones found by Nakayama [38]. These two oriented 
matroids can also be realized using Grobner basis. 

Theorem 4.1 There are precisely 1, 1 and 2 irrational realizable oriented matroids in OM(3, 9), OM(6, 9) 
and OM(4, 8) respectively. 

From these results, we obtain the combinatorial types of point configurations (Theorem 1 1.2| by gen- 
erating relabeling classes of acyclic realizable oriented matroids. Matroid polytopes (i.e., acyclic oriented 
matroids with all elements extreme points. For details, see [5].) are extracted from them. Then we 
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compute the face lattices of the matroid polytopes and decide whether they occur from some realizable 
matroid polytopes or not in order to obtain combinatorial types of polytopes (Theorem II. 2[) . All face 
lattices of matroid polytopes in OM(4, 8),OM(3, 9) and OM(6,9) turned out to be realizable as those of 
convex polytopes. Table [7] and Table |S] summarize the results. 





acyclic OMs (realizable/non-realizable) 


matroid polytopes (realizable/non- realizable) 


OM(4, 8) 


10, 775, 236 (10, 559, 305/215, 931) 


250,601 (238,399/12,202) 


OM(3,9) 


15, 296, 266 (15, 287, 993/8, 273) 


1 (1/0) 


OM(6,9) 


105, 183, 749 (105, 128, 749/55, 000) 


41, 030, 709 (41, 008, 968/21, 741) 



Table 7: Numbers of relabeling classes of acyclic OMs and matroid polytopes 





total (realizable/non-realizable) 


OM(4, 8) 


257 (257/0) 


OM(3,9) 


1 (1/0) 


OM(6,9) 


47,923 (47,923/0) 



Table 8: Numbers of non-isomorphic face lattices of matroid polytopes 

The number of non-isomorphic face lattices of 3-polytopes with 8 vertices coincides with the number in 
[3"ffl p. 424, Table 2]. We observe that one can associate a point configuration with rational coordinates 
to the combinatorial type of every 5-polytope with 9 vertices and thus obtain the following theorem. 

Theorem 4.2 The combinatorial type of every 5-polytope with 9 vertices can be realized by a rational 
polytope. 



5 Concluding Remarks 

In this paper, we complete the realizability classification of OM(4, 8), OM(3, 9) and OM(6, 9) by develop- 
ing new techniques to search for a realization of a given oriented matroid. Surprisingly, the biquadratic 
final polynomial method [5], which is based on a linear programming relaxation, can detect all non- 
realizable oriented matroids in these classes. In addition, one can also find all non-realizable uniform 
oriented matroids in OM(3, 10) and OM(3, 11) by this method [31 0]. A known minimal non-realizable 
oriented matroid which cannot be determined to be non-realizable by the method is in OM(3, 14) [4"2"j . 
It may be of interest to find a minimal example with such property. 

Our classification almost reaches the limit of today's computational environments. However, we can 
deal with larger instances if the number of instances is small. Actually, it was successfully applied to 
OM(4, 9) and OM(5,9) in order to find the hyperplane arrangements maximizing the average diame- 
ters [T7] and PLCP-orientations on 4-cube [TH]. Our classification results are available at 
http : //www-imai . is . s .u-tokyo . ac . jp/~hmiyata/oriented_matroids/ 
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